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Abstract 

A new model method for describing of the electrostatic screening in two-component systems 
(electron-ion plasmas, dusty plasmas, electrolytes, etc) is developed. The method is applicable 
to the systems of higher non-ideality degree. The expressions for all the screening parameters 
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characteristic for multi-component systems, are obtained. All these parameters are presented in a 
simple analytic form suitable for operative laboratory usage, especially for theoretical interpretation 
of experimental data. 
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I. INTRODUCTION 

In the previous paper here Part 1, the aims of this research were already described, 
as well as the stimuli for its starting. Accordingly to these aims the new model method of 
describing of the electrostatic screening in electron-ion plasmas and other two-component 
systems (e.g. dusty plasmas and some electrolytes) which relies on the basic model (al) 
- (a3) is presented in this part. This method is free of non-physical properties of Debye- 
Hiickel's (DH) method and posses the positive features (bl) and (b2) described in Section 
1 of Part 1, together with the basic model. 

The material presented in this paper is distributed in the four following Sections and four 
Appendixes. Section [III contains: the screening model; the critical analysis of DH method in 
the case of two-component system; stating the tasks precisely. In Section UTTl and Section HVl 
the developed method and obtained solutions for two-component systems are presented. 
The Section |V] contains obtained results and discussion. 

II. THEORY ASSUMPTIONS 
A. Screening model 

A stationary homogeneous two-component system is taken here as the initial model 
of some real physical objects. We will assume that Sm is constituted by a mix of two gases: 
one of positive charged ions (of only one kind), and other of electrons. It is assumed that 
these gases there are in the equilibrium states with temperatures Tj and Tg > Tj, and mean 
local particle density Ne and Ni. All the particles are treated as point objects with the 
charge ZgC in the case of electron, and Z^e in the case of ion, where = —1, Zi = 1,2,..., 
and e is the modulus of the electron charge. Let us note that in this paper the electron 
charge will be also denoted by — e. It is understood that the parameters Ze,i and A^e.i satisfy 
the local quasi-neutrality condition 

Zie-N,-e-Ne = 0, (2.1) 

as well as that Ng and Tg allow the non-relativistic treatment of the electron component. 

In accordance with the properties (al) and (a2), the screening of a charged particles in 
the system Sin will be modeled in the corresponding accessory systems each of which differs 
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from Sin in that, besides the two described components, it also contains a fixed probe particle 
with charge ZpC in the origin of the chosen reference frame (point O). Here we will study 
two cases: the ion case (i), when Zp = Zi, and the electron case (e), when Zp = Z^. = —1, 
when the probe particle represents one of the particles of the system Sin- In accordance with 
this, we will denote here the corresponding accessory system with Sa^ or s!f\ This system 
will be characterized by: the ion and electron densities nf'^\r) and ne'^\r), the mean local 
charge density 

p(i.-)(r) = Z,e ■ nf'^^ir) - e ■ n(*'^)(r), (2.2) 

and the mean electrostatic potential $(*'^)(r), where r = |r| and f^is the radius- vector of the 
observed point. It is assumed the satisfying of the boundary conditions 

lim nf'''\r) = Ni, Mm n'^^'^\r) = N^, (2.3) 
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and the condition of neutrality of the systems Sp''^'' as wholeness 

oo 

Z„e + /p<«'M.W*^0. (2.4) 



Then, we will take into account that the <|>(*'^)(r) and have to satisfy Poisson's 

equation 

V2$(^.-) = _47r [Zi^ee ■ 5{r) + p^'''\r)] , (2.5) 

n 

where (5(r) is three-dimensional delta function [2]- From the same reason as in Part 1, this 
equation applies in the whole region r > 0. It is assumed the satisfying the boundary 
conditions 

lim $(^'^)(r) = 0, (2.6) 



|(/?(^'")| < oo, = lim[$(^''=)(r) - Zi^e/r]. (2.7) 

Since y?'^*''^) is the mean electrostatic potential in the point O, the quantity 

U^''^'^ = Zi^.e ■ (2.8) 

is the potential energy f/(*'^) of the probe particle. In an usual way U^^^ and Z7'-^^ are treated 
as approaches to the mean potential energies of the ion and electron in the system Sin- 



In accordance with the properties (a2) and (a3) the conditions of thermodynamical equi- 
hbrium of the ion component in the case {i), as well as the electron component in the case 
(e), will be taken in the form 

^^^,e (nff{r),T,,,^ + Z,,ee ■ ^^'^'\r) = fa,, (iVi,e,T,,e) (2.9) 

where fii{nf\r),Ti) and ^e{n^e\r),Te) are the chemical potentials of the ideal ion and elec- 
tron gases, which can depend of the corresponding particle spins, considered on the distance 
r from the point O. On the base of the considerations from Part 1 one should keep in mind 
that the equations (12. 9p are applicable only in the regions 

1/3 



r > r. 



■,i,ei 



(2.10) 



where Vg-i and rs;e are the corresponding Wigner-Seitz's radii. In the used procedure the 
equation (12. 9p is taken in the linearized form 



-N- = 



but under the condition 



dn 



n=Ni, 



(2.11) 



< 1. 



(2.12) 



It is important that the conditions (I2.10p and (12.120 are compatible in all considered cases. 

Since we take the single-component systems considered in Part 1, as a boundary case of 
two-component systems (when it is spread one of their components), we will require that 
the ion density nf^ in the case (i) and the electron density rif^ in the case (e) satisfy the 
equation 



Anr'^dr = 1. 



(2.13) 



which is analogous of the equation (21) from Part 1. From Eqs. ( 12. 4p and ( 12.13^ it follows 
that the electron density ni*^ in the case (i) and ion density nf''' in the case (e) have to satisfy 
another equation 



■ Anr'^dr = 0. 



(2.14) 



Let us emphasize that the relations fl2.13p and (12.141) can be obtained on the base of the 
interpretation of the systems Sa''^^ which is given in Appendix [Dl In further considerations is 
used the fact that simultaneous satisfying of the conditions ( I2.13P and fl2.14p automatically 
provide the satisfying of neutrality condition fl2.4p . 

B. The critical analysis of DH method 

The procedure of obtaining of DH solutions is described in Appendix El The figure [H 
shows the behavior of the particle densities n^^.-{r) and n^^.^{r), and charge density p'j^{r). 
This figure illustrates that the behavior of DH densities of the free particles with the same 
charge as the probe particle is qualitatively same as in the case which is considered in the 
Part 1. Consequently, the procedure of the elimination of the non-physical properties of 
n^^^(r) in the case (i), and n^^.f.{r) in the case (e) will be similar to the procedure which is 
described in Part 1. 

The main disadvantage of DH method consists in the monotonous increasing of DH 
densities of the free particles, which the charge is opposite to the charge of the probe particle, 
with the decreasing of r in the whole region r < oo. This fact is illustrated by the behavior 
of Ti~Q.g^[r) in Fig. [1] Because of such a behavior DH solutions principally can not satisfy 
the conditions (12.140 . Namely, from flASp it follows that in DH case the left side of those 
conditions is not equal to 0, but it is proportional to Zi in the case (i), and (1/^j) in the case 
(e). The consequence of this fact is the principal impossibility to treat the probe particle as 
a represent of a particle in the system Sin- For an example, in the case of completely classical 
electron-ion plasma with Zi = 1 and Tj = from the non-satisfying of the conditions (I2.14p 
it follows that the mean number of electrons per ion should be 3/2, instead of 1. 

The described disadvantage is a consequence of two facts: that in DH method the first 
step is determining of electrostatic potential $^'^'*(r) in the whole space, and the equations 
(12. lip are used together with the equations flAip from Appendix |X1 Because of that within 
DH method the electron and ion components are "smeared" in the space simultaneously and 
independently in both (i) and (e) cases. 
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r [arbitrary units] 



FIG. 1: The reduced DH densities n^^.-{r)/Ni, n^^.^{r)/N(. and {r) / {—eNf.) in the case Zi = 1, 
Te = Ti and K£)rs-i = 1, where kd is Debye's screening constant given by (|A2p . Shadowing 

(i) 

emphasize the non-physical deviation of njj.^ from the asymptotic value A'^e- 

C. What one should do in order to avoid non-physicality of DH method? 

In accordance with above mentioned, our main task is finding of such a procedure within 
the basic model (see Section 1 in Part 1), which would be alternative one to DH procedure. It 
assumes that sought procedure has to provide the possibility of determining of the solutions 
nf\r) and ne\r) in the same way as it was described in Part 1, and the solutions ne\r) 
and n^e\r) without using the equations (lAip from Appendix |Xl Also, the satisfying of the 
conditions fl2.13l) and fl2.14p is assumed. 
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III. THE PRESENTED METHOD: THE ELECTRON AND ION DENSITIES 



A. The case (i) 

The solution nf\r). In order to solve our task we start from such a picture of the 
system Sa'' where electrons are treated in the electrostatic field of the probe particle and 
all ions (treated in a classic way) which are distributed in discrete points, as it is illustrated 
by Fig. [21 On the base of this picture, the procedure of the expressing of n^e\r) in the 
region r > r^.j trough nf\r) is developed in the first part of Appendix [Bl The results of 
this procedure is the relation flB3p . which can be presented in the form 

n^^\r)=N^-{l-a) + a-Zinf\r), r..^ < r < oo, (3.1) 

where the parameter a will be determined in further text. This relation, together with 
Eq. (12.21) . makes possible to represent the charge density p^^\r) in the form 

p^'\r) = Z,e ■ {1 - a) ■ [nf\r) - Ni], r,., < r < oo, (3.2) 

and to use for determination of nf\r) the procedure which is described in details in Sections 3 
and 6 of Part 1. As first, by means, Eqs. (12. lip and (13.21) it is obtained the equation of 
Volterra's type, namely 



oo 













)(r') 










\r 


r'J 



pi^{r) = kI / p«(r') - - -\r'Hr', (3.3) 



where 

Ki = = KO;i ■ (1 - a) / , KO;i 



(3.4) 



^^i^|^N,_ 

One can see that in the two-component case the ion screening constant Hi and the corre- 
sponding characteristic length r^jj depend on the parameter (1 — a)^. 

Accordingly to Part 1, the solution of the equation (13.31) in the general case can be taken 
in the form: p^^\r) = Ci-Zie-exp(—Kir)/r. From here and Eq. (13. 3p it follows that ion density 
nf^(r) in the region Vs-i < r < oo is given by the relation: nf\r) = Ng — Ci ■ exp(— Kjr)/r. 
Than, taking Cj = r^-i ■ exp(Kjro;i) ■ A^e and applying to this relation the extrapolation 
procedure from Part 1, we obtain the sought expression for n'f\r) in the whole space. 
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FIG. 2: The schematic picture of ion-electron plasma. Shadowed areas represent the regions of r 
without of probe particle's and all ion's self-shells. With O and O' the places of the probe particle 
in the case (z) and one of ions are denoted, and with the dotted lines probe particle and ions self 
spheres are shown. 

namely 



AT AT / \ exp(— Kjr) 

n\ '{r) = { r (3.5) 

0, r < ro;, 



where the new characteristics length r^-i has to be determined from the condition (12.131) . 
Using the results of Part 1 we can present the radius r^-i in two equivalent forms 

ro;i = 7s{xi) ■ Ts-i, ro-i = -f^{xi) ■ r^-i, (3.6) 



(3.7) 



where the coefficients 7s(a;) and 7K(a;) are given by Eq. (28) from Part 1, and Vg-^i - by 
Eq. 

The solution ni\r). In order to determine rii\r) in the probe particle self-sphere 
(0 < r < r^.j), we will take into account that in the system Sin the mean number of electrons 
in the sphere with volume which is centered at some ion (the ion self-sphere), is 

larger than the mean number of electrons in every fixed sphere with the same volume, i.e. 
Ne ■ iX/Ni) = Zi, because of the additional electrons whose presence is caused by presence 
of other ions. In accordance with this we will find ni* ■* (r) in the probe particle self-sphere in 
the system Sa^ in the form 

nW(r) = ne;s{r) + ne-ion{r), < r < r^.^, (3.8) 
where the member ne-s{r) satisfies the condition 

j ne;s(r) ■ Airr^dr = Zi, (3.9) 



considering that this member describes the distribution of Zi electrons which there are in the 
probe particle self-sphere independently of the presence of ions, while the member ne;ion{f) 
describes the distribution of the mentioned additional electrons. Let us emphasize that such 
a treatment of the difference [ne\r) — ne-s{r)] = rie-ionir) is caused by the procedure of the 
obtaining the member ne;s(r) which is described in the second part of Appendix fiBl) . 

In order to establish the connections between the members in Eqs. (13.11) and (13. 8p . we 
will assume that the form of the equation (13.11) . which transforms at r = oo to equality 
Ne = Ne{l — a) + NeO-i reflccts the existing in the system Sin of the spatial correlation 
between electrons and ions. It means exceptionally following: in every fixed volume V in 
average there are V ■ Nf.{l — a) electrons independently of presence of ions, and V ■ N^a 
electrons whose presence is caused by the presence of ions. The fact that the presence of ions 
in some volume causes the presence of additional electrons itself has already been used above, 
while here it is assumed that just the parameter a in Eq. (13. ip represents the quantitative 
characteristic of the mentioned electron- ion correlation in the system Sa\ Because of that 
the first of the mentioned connections is given by the relation 

ne;ion{r) = aZi ■ nf\r), < r < r^-i, (3.10) 
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which corresponds to the the way of obtaining of the ion density, and provides the satisfaction 
of three conditions: that the member ne-ion{r) in (13. 8p has to be close to aZi ■ n[^\r) at least 
in one part of the region r < r^.j; in whole this region the behavior of Ue-ioni''^) has to reflects 
the behavior of nf\r)] the ratio of the mean numbers of additional electrons and ions in 
the probe particle self-sphere has to be equal to correlation coefficient, i.e. a under adopted 
conditions. Another of the mentioned connections is the condition 

ne■,s{rs■,^) = N,{l~a), (3.11) 

which provides needed continuality of ni\r) at r = r^jj. The physical sense of this condition 
is discussed in the further text. 

The procedure of determination of the member ne;s(r) is described in the second part of 
Appendix [Bl The first results of this procedure are the equation flB4l) and the expression 
( IBSp . which provide the obtaining of the equation for direct determination of ne-s{r), namely 



^0;e o /" , / 1 1 

Att r 





4vr Tst J \r' Tst 



where the screening parameter kq-^. is given by 



(3.12) 



1 



The next results of the used procedure are the relations (IBlOp and (IBllI) which determined 
the member UfAr) in the form 



r ' \ rs-ij r ' \ r^.. 



a ■ — '- exp I —Xs — ] + ■ — - exp I Xg — 



0<r<r,.i, (3.14) 



where the coefficients a and b have to satisfy the condition fIBlSp . which provides that this 



relation really represents a solution of the equation (I3.12p . as well as the condition (]B16p . 
which provides the satisfying of the condition (13.111) . Under these conditions a and b are 
obtained in the form 



1 — a — exp(xs) 1 — a — exp(— a;^ 



(3.15) 



exp(-Xs) - exp(xs) ' exp(-Xs) - exp{xs) 

= l^0;ers;i- (3.16) 
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Finally, the last free parameter, i.e. the correlation coefficient a, is determined from the 
condition (13. 9p and is given by the expression 

2^,3 

« = l-7^ ^ 1 ^ r^- (3-17) 

(1 + Xs) exp (-Xs) - (1 - exp (xj 

From here it follows that a satisfies the conditions 

< a < 1, < < cxd; lim a = 0; lim a = 1; (3.18) 

a;s— >0 Xs— >oo 

which make possible the treatment of a as the correlation coefficient in the whole region 
< < cxD. Finally, by means of Eqs. (13. ip . (13. 8p . (I3.10p and (13. lip the electron density 
ne*^(r) can be presented in the form 

A^e(l - a) + OiZi ■ nf^(r), r^.j < r < cx), 
^e;s('^) + ttZj ■ nf'[r\ < r < 

where nf''(r) and r;,e;s(r) are given by Eqs. (13.50 and ( 13.14p - (13.16p . which provide ne*^(r) 
satisfies the condition Eq. (I2.14p . while Oj and 6i are given by (I3.15P and a- by (I3.17p . 

It can be shown that ne;s(r), given by (I3.14p - (l3.17p . in the region r < Tg monotonically 
increases with decreasing of r and satisfies the equality 

^^^^^ =0, 0<x, <oo. (3.20) 

Apart of that, these facts provide smoothness of the electron density rSe\r) in the point 
r = Tg-i and give the clear physical sense to the condition (13. lip within the developed 
method. Namely, we have that just the value ne-^s{fs;i) in this condition determines the value 
of the effective density of the non-correlated part of the electron component, i.e. Ne ■ (1 — a), 
inside and outside of the probe particle self-sphere. 

Finally, the property (I3.20p of r2e;s(?") causes such a behavior of the difference [ne-si.'^) — 
^e;s('"s;j)]5 whlch is illustratcd by Fig. [31 This figure shows that the deviation of that difference 
from zero can be practically neglected within the layer 0.75rs;j < r < r^.j for any x^ > 0. 
This provides the extrapolation of the rLi\r) from the region r > r^^j in the part of region 
r < Ts-^i which makes about 60% of the probe particle self-sphere can be performed in the 
same way as in Part 1, which means that the deviation of [ne;s(r) — ne;s(rs;j)] from zero can 
be neglected. It is clear that such a behavior of Ue-sci''^)-, which could not be expected in 
advance, additionally justifies the applied procedure. 
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FIG. 3: The behavior of the ratio [ne;<i(r) — ne-s{fs\i)]/Ne^ for 0.1 < Xg = KQ-e'>'s;i < 20. 
B. The case (e) 

The solution n^e\r). In this case treatment of the light (electron) and hard (ion) 
components appears as the main problem. The way of its solving and the procedure of the 
expressing n^f\r) through n^e\r) in the region rs-e < r < oo are described in first part of 
Appendix O The results of this procedure is the relation fICip which can be presented in 
the form 

n'f\r)=Ni{l-a) + ^-ni^\r), r,;e < r < oo. (3.21) 

This relation, together with Eq. 02.21) . makes possible to express the charge density p^^\r) in 
the form: p^'^^r) = —e ■ [nit\r) — N^] ■ (1 — a), and repeat word-for-word the procedure from 
Section IIII A[ As result we obtain the needed expression for the electron density ni^^ (r) in 
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the whole space, namely 

{AT AT / \ exp(— Kgr) 

Ne - N^ro-e ■ exp(Kero;e) , ro-e < r < oo, 

r (3.22) 

0, < r < ro;e, 

= — = fi:o;e ■(!-«) ^ (3.23) 

where kq-s is defined by (13.131) . and the electron screening constant Kg and the corresponding 
characteristic length ri^.^ depend on the parameter (1 — a)^/^. The radius ro;e, similarly to 
the case (i), can be presented in two equivalent forms 

rO;e = 1s{Xe) ■ rs;e, TQ-e = ^K^Xe) " T^-e (3.24) 



L~ If 



(3.25) 



where 7s(a;) and 7k(x) are given by relations (28) from Part 1, and Vs-e - by Eq. (I2.10p . 

The solution n'f^ (r) . Repeating the procedure from Section IIII Al we will find the ion 
density n['^\r) inside the probe particle self-sphere, < r < r^-e, in the form 



n'f\r) = fii-sir) + rii-eiir), rii-eiir) = ^ ■ n^eK^) 



a 



(3.26) 



where the member ni;s(r) satisfies the conditions 



rii-sir) ■ 47rr dr = 



(3.27) 



ni.s{rs;e) = Ni{l-a), (3.28) 

which play similar role as the conditions (13. 9p and (13. lip in the case (i). 

The way of determination of rii-sir) is described in second part od Appendix [Cl and as 
the result the relations flC2l) and flC3P are obtained. After determination of the coefficients 
in the superposition (lC2p from the conditions (I3.27P and (13.280 . the member nj;<j(r) can be 
presented in the form 



rse I ^ \ 1 '^'s-e / f 

a ■ — '- exp —Xs + ■ — - exp Xs 

r V r,.p / r V r,.^ 



, < r < Ts-e, 



(3.29) 
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where a and b are given by Eq. fl3.15l and Xs - by Eq. (13.161) . Finally, from Eqs. (I3.2ip . (I3.26P 
and (I3.29P it follows the expression for n['^\r) in the hole space, namely 



Ni{l -«) + — ■ nf>{r), rs-e < r < oo. 



(3.30) 



where ni'^\r) and rii-sif) are given by Eqs. (I3.22p - (l3.25p and (I3.29p . which provide that 
n[^\r) satisfies the condition (12.140 . and a is given by Eq. (I3.17p . 

It can be shown that the member ni-s{f) monotonously increases in the region r < Vg-e 
with the decreasing of r and satisfies the equality 

dni-s{r) 



dr 



0, 



< X, < CXD, 



(3.31) 



which provides smoothness of rLf\r) in the point r = r^.g. Consequently, it guarantees that 
the ion density rtf\r) in the case (e) has the similar properties as the electron density n^e\r) 
in the case (i). 

IV. THE SOLUTIONS p(^ ^)(r) AND $(^ '=)(r), AND THE PROBE PARTICLES PO- 
TENTIAL ENERGIES ?7(*'^) 



By means Eqs. (Q, O, fl3J9|) . (13:221) and fl330D the charge densities p(*'^)(r) can be 
presented in the form 



(4.1) 



< r < cxo, 



(4.2) 



< r < r. 



(4.3) 



where n^le\r) and ^^e*^^s(^) given by Eqs. (13. 5p and (I3.22p . and Eqs. (I3.14p and (13.290 . 
respectively. One can see that the expression (14. 2 p has the same structure as the correspond- 
ing expression from Part 1 and becomes identical to it for a = 0, while (14. 3 p describes the 
contribution of the quantity \n!'e,i-!si^) ~ — a)], which characterizes just two-component 

systems. 



14 



In accordance with the structure of the expressions for p(*''^)(r) we will find the electrostatic 
potentials $(*'^)(r), as well as the potentials yj^*'^) defined by Eqs. fl2.7lj . in the form 



(4.4) 



where the first and second members describe the contributions of the members pa''^\r) 
and pa-t\r) in Eq. (14.11) . It can be shown that by means Eqs. (14. 2 p and (14. 3p . as well as 
Eqs. (51)-(53) from Part 1, the members ^a'^\r) and ^a-s\r) can be presented in the form 



Zi,e-{l-a) xKe) ■exp(-fi:,,er), 

(^2''V 1 

1 Ji \- _ 

Ziee-{l-a) 2 \r. 



ro-i,e <r < oo, 
<r < ro;i,e, 



(4.5) 



a + — h 3 • — 



) + 6exp(xs-^) + (a - b)xs 



Zi^f,e 1 — a 



< r < r^.j 



i,ei 



- 1 



(4.6) 



where Xj and Xe are defined by Eqs. (13. 7p and (I3.25p . Xg - by Eq. (I3.16p . the function x{x) 
- by Eq. (35) from Part 1. The quantities ipa''^^ and (falt^ are obtained by means Eqs. (14. 2 p 
and (14. 3p . as well as Eq. (52) from Part 1, and presented in the form 



y^i^"^^ = -^.ee ■ (1 - «) ■ 



3ro;i,e / , 2 

rO;i,e + - 
K 



(4.7) 



ra:s 



2r • 



ail — e 



+ 6(e^» - 1) 



-il-a) 



Xa 



(4.8) 



where the coefficients a and h are given by Eqs. (13.150 . and the radii ro;i,e - by Eqs. (13. 6p . 
<K7\ . (13:2^1) and fl3:25|) . as well as Eq. (28) from Part 1. 

In accordance with Eqs. (12. 8p . (14. 4p . (14. 7p and (14. 8 p it is appropriate to find the potential 
energies f/^*'*^) of the probe particles in the cases (i) and (e) in the form 



(4.9) 



where the quantities Lfa''^'^ and ifa-f are given by Eqs. (14. 7p . (14. 8p . It can be shown that, 
because of the structure of Eq. (14. 7p . the member Ua''^^ can be taken as: Ua''^^ = {l — a)-U, 
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with U given by one of two equivalent expressions (32) and (33) from Part 1, where the 
parameters Z, x, and are replaced by Zj^e, Xi^e, rK-i,e and rg-i^e- The figure (jl]) illustrates 
the behavior of the ratio of the probe particle potential energy, determined by means to 
Eqs. fl4.7p . (14 .Sp and (14.91) . and the corresponding DH energy, defined by Eq. flA6p . in the 
case of the classical plasma with Zi = 1 and Tg = Tj. This figure shows that in two- 
component case the DH region, i.e. the region where this ratio is close to unity, principally 
does not exist. 



V. RESULTS AND DISCUSSIONS 



The expressions (13. 50 - (13. 190 and (13. 220 - 04. 5p show tha.t:the obtained solutions satisfy all 
conditions from Section are free of all non-physical properties of DH solutions, and pos- 
sesses the positive properties (hi) and (h2), noted in Section 1 of Part 1. The behavior of 
electron, ion and the charge density is illustrated by Fig. [5l The way of their obtaining 
provides the possibility of their application, independently of the treatment (classical or 
quantum mechanical) of electron and ion components of the considered systems. Apart of 
that, the structure of these expressions makes possible their application not only to electron- 
ion plasmas, but also to ion- ion two-component systems (e.g. some electrolytes). 

Since Eqs. (I3.14p . (13.190 . (I3.29P and (I3.3O0 show that the solutions nl^f\r) are singular 
in the point r = 0, it is useful to note that the existence of singularities in model solutions 
is fully acceptable, if it has not other non-physical consequences. Such solutions are well 
known in physics: it is enough to mention, for example, Thomas-Fermi's models of electron 
shells of heavy atoms (Q, 4|; see also [sl), which use in plasma research till now (see e.g. 



In the case of two-component system are obtained the parameters ro;i,e, Isi^i^e) and 
luixi^e), given by Eqs. (13. 6p . (13. 7p . (I3.24p and (I3.25p . analogous to that ones from Part 1. In 
accordance with Eq. (28) from Part 1 the parameters r^.i and tq-^ satisfy the conditions 

< ro;i,e < rs-i,e, < Xi-e < OO] lim ro;i,e = 0; lim ro;i,e = Vs-i^e] (5.1) 

where r^.j e is given by Eq. (I2.1O0 . and Xj and Xe - by Eqs. (13.70 and (I3.25p . These conditions 
make possible the treatment of tq-i and ro;e as the radii of the spheres centered on the 
probe particles, which are classically forbidden for the free charged particles from their 
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0;i s;i 



FIG. 4: The behavior of the ratio as a function of the parameter Ko-trg-i in the case of 

completely classical plasma with Zi = 1 and T, = Tg, where KQ-i is defined by Eq. (|3.4p . 

neighborhoods, and 7s(a;j^e) and 7K(xi^e) - some kind of non-ideality parameters (see Part 
1). From it follows also that Eqs. ([33]), fl3J9D . fl3:22|) and (gS]) for electron and ion 
densities are applicable to the two- component systems with any non-ideality degree. 

Also, in this paper is obtained the quantity a{xs), defined by Eqs. (13.161) and (13.171) . 
which has the sense of the coefficient of electron-ion correlation. Let us note that two simple 
approximative expressions for a{xs), which serve very well in wide region of Xg, are given in 
Part 3. 

One of the most important results of this papers is establishing of the fact that in two- 
component plasmas ion and electron components have to be described exceptionally by 
means of screening constants Ki and Kg, and the corresponding screening radii r^^-i and r^jg 
which are introduced in this paper. This means that Debye-Huckel's screening constant kjj 
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FIG. 5: The behavior of reduced densities n'f\r)/Ni, ni\r)/Ne and p^'^\r) / {—eNf.) in the case 
= 1, Tg = Tj and KQ-iVg-i = 1, where kq-j is defined by (|3.4p . 



and t/ie radius ro do not appear in the theory and, consequently, principally do not have the 
physical sense. It is confirmed by Fig. Hlwhere the probe particle potential energies obtained 
here are compared with the corresponding DH values. Finally, comparing of the expressions 
f l4.4p . fl4.5p and (14.61) for the potentials $(*'^)(r) with the expression flA3p for DH solutions 
$^'^''(r) shows that in two-component case the principal difference between these solutions 
there is not only inside the probe particle self-spheres, but also in the rest of space. Namely, 
out of these self-spheres $^'^''(r) ~ exp(— /t£i)r, where kd is given by Eq. flA2p . while the 
solution $(*''=)(r) ~ exp(—Kj^e)?% where Ki and Kg are given by Eq. (13. 4p . (13.130 and (13.230 . 



This fact justifies the usage in 
of DH constants. 



10| the constants which are close to and Kg instead 
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APPENDIX A: DH SOLUTIONS 

In DH method it is assumed that for determination of rie^ and nf\ apart of Eqs. fl2.11l) . 
it can be applied the equations 



g/ie,t(n,Te 

dn 



(Al) 



which represent the linearized form of the equations: ^e,i (^^e f\''^)y'^e,i^ + Z^^iC ■ <|)(*''^)(r) = 
fJ^e,i {Ne^i,Te^i). Within DH method equations fl2.1ip and flAip apply together in order to 
express the charge densities in the Poison's equations (12. 5p through the electrostatic potential 
in the whole region < r < oo. That way, in accordance with Eqs. (I2.1l) - (l2.3p and (12.51) one 
obtains Helmholtz's equation 

V'^(^'^\r) = ^D^- = {4, + <e)K (A2) 

where kd is DH screening constant, and Ko;j and Ko;e are the partial screening constants 
defined by Eqs. ([32D and fl37[3D . 

DH electrostatic potential $^'^^(r) represent the solutions of equation flA2p . which is 
obtained in the whole space under the boundary conditions (12.61) and (12.71) . Then, by means 
of Eqs. (12. lip and flAip DH solutions for the charge and particle densities are obtained. All 
these solutions are given by the expressions 

\n = exp(-Kz)r), p'^ '{r) = — , (A3) 



= + „g) (,) = ^v. + ^^2^. (A5) 
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where the screening constant is given by (lA2p . From Eqs. (12. 7p . (12.80 and (1A3P it follows 
that DH potential v^^''^'' = — ^i,ee ■ and consequently, DH potential energy U^'^^ is given 
by 

U^'^'> = -Z,,,e-^n = -^-^^, (A6) 

rD 

where and are defined by Eq. (lA2p . 



APPENDIX B: THE SYSTEM si*^ 

The region r^.j < r < oo. We will start from the fact that the different electron-ion and 



dusty plasmas can be successful 



particles (see e.g. 



12 



13 



14 



l y d escribed in the approximation of fixed heavy charged 



151]). In accordance with this, we will treat the electronic 
component of the system Sa^ considering all ions as immobile with respect to electrons, 
keeping in mind that in such a case the difference between the probe particle and ions does 
not exist regarding the electrons. In Fig. [2] the several of ions self-spheres (sphere with the 
radius Vg-i, centered on ions) in the neighborhood of the probe particle (the point O) are 
schematically shown. This figure should make more apparent the fact that the behavior of 
the electron inside the self-sphere of an ion in the point O', far from the point O, has to 
be almost completely caused by its interaction with few ions which are close to the point 
O'. Because of that, in further considerations all ions are treated equally with the probe 
particle. 

Here, we will denote by (*'*) the possible ion configuration in the system Sa\ and by 
ne'*\r) and n'f'*\r)- the corresponding electron and ion densities. Also, we will denote by 
Ue-'olt average value of Ue'*'^ (r) within the part of space consisting of the outing of probe 
particle's and all ion's self-spheres (see shadowed area in figure [2]). The condition of the 
thermodynamical equilibrium of the electron component in the case of a configuration (*•*) 
has the form 

/ie (ni^^*Vl,Te) + (-e) ■ $(^^*)(r) = /i, ^^^^^(rl*), T^) + (-e) ■ ^^^*^{f^:), (Bl) 
where (r) is the corresponding total electrostatic potential, and f^^'*'' is any point where 



n. 



e'*\r) = n^e-^*l^. Finally, taking the chemical potentials fie{ne'*\r),Te) in (IBll) as series of 
the differences [rie'*^ (r) — Ne] along with keeping only first two members, we obtain ne*'*'' (r) 
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in the form 



— e 



In contrast to (12. lip , the equation (lB2p can be apphed practically in the whole space, which 
is allowed by the behavior of the electron component in the presences of the positive charged 
particles. 

In accordance with the described picture, the sought electron and ion densities n'e\r) and 
rif'{r) are treated here as a result of averaging of all densities rti'*\T) and ^f'*^ Because of 
the structure of Eq. (lB2p we will have that: ne\r) = Ng-st + 'ne-d{i^), where N^-st and ne-dir) 
are the mean values of the first and second members in Eq. flB2l) . Taking the member ne;d(r) 
as series of the differences [nf \r) — Ni] along with keeping only first two members, we obtain 
that: ne\r) = Ne-st + ne-dir = oo) + K ■ [nf\r) — Ni], where accordingly to the boundary 
condition (12. 3p the member ne-d{r = oo) has to satisfy the equality ne-d{r = oo) = Ne — Ne^st- 
Finally, taking that K = a ■ Zi we obtain the mean electron density ni* ■* (r) in the form 



n 



Ne + a-Zi[nf>{r)-Ni]. (B3) 



This relation shows that ne\r) in the region r > Vg-i can be expressed by means only one 
unknown parameter a. 

The region < r < r^.j: the basic equation. Let (^'*) denotes such an ion configura- 
tion in which the probe particle self-sphere is left free of ions, and ni^'*^ (r) - the corresponding 
electron density. The member ne;s(r) in (13.80 will be identified with the result of averaging 
of densities ni^'*\r) over all configurations ('''*) in the region < r < rg-i- 

The corresponding condition of the thermodynamical equilibrium in the case of a con- 
figuration ('''*) obtains from (IBip replacing the index (i,*) by (s,*) and fixing the point 
fjH'*^ = fst, where Tgt < r^^j. Then, repeating the procedure of the obtaining of Eq. flB2p . 
we obtain the equation 

rie-Ar) - rie-Arst) = q^j^j^^ " [^s{rst) - $.(r)] . (B4) 

where the potentials $s(r) and ^si'f'st) in the considered region are given by Eq. (53) from 
Part 1 with the charge density (— e) ■ rie-sif^)- Since all information about the outing of 
the probe particle self-sphere {vg-i < r < oo) contained in the same constant members in 
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expressions for $s(r) and $s(rsf), the difference [^si'f'st) — ^s{f)\ Eq. (lB4p is given by 









y (-e) ■ n^-s{r 










rst 



+47r y (-e) ■ - r'^cir'. 



(B5) 



From Eqs. (lB4p and (]B5|) it follows the equation fl3.12p for direct determination of the 
member r;,e;s(r). In order to find the solution of the equation (13 ■121) . we will introduce the 
function Sir) given by relation 



n , r 



Air r r 

Consequently, Eq. fl3.12jl . after the multiplication by r, transforms to the equation 
S{r) — Kp.g / S{r') (r — r') dr' \_^''^ / ('^ ~ ''^O c^'^' 



(B6) 



47r 



r 



^(r^i) - Ko-e / S{r') {rst - r') dr' - 



' ^0;e 

47r 



(B7) 



Applying the operator ^ to left and right sides of IB 71 we obtain the equation 



d^_ Zj^ 

which can be presented in the form 
d^ 



0, 



(B8) 



dr'^ 



Sir) 



ZiK^.^ 

An 



Sir) + 



ZiKq.^ 
At: 



(B9) 



From Eqs. (]B6P and (]B9P it follows that in the general case ns{r) is given by the relation 
A ■ exp {-KQ-er) + 5 ■ exp {nQ-er) 



risir) 

r 

Taking the coefficients A and B in the form 
A = rs-i ■ Ne- a, B = r^-i ■ Ne ■ b, 



(BIO) 



(Bll) 
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we will present Eq. (13.121) in the form 
a ■ exp (— Ko;e'") + & ■ exp {^Q-^fS) 



+Ko.g y [a ■ exp {-KQ-er') + & ■ exp {n^-er')] ■ (r - r') dr' 



= — < a ■ exp {-Ko-erst) + b ■ exp (Ko.erst) - 

^st I 47rrs.jiVe 

rat 

+Ko.g J [a ■ exp {-Ko-er') + h ■ exp {i^Q-er')] ■ {rst - r') dr' 



(B12) 



From here it follows the equation 
(6 — a) KQ-eT + I a + 6 — 



r 



{h - a) Ko-eTst + [a + h- 



(B13) 



which can be presented in the form 

a + 6-^^ = ^fa + 6-^^), (B14) 



where in the general case r ^ r^^. This means that the coefficients a and b have to satisfy 
the condition 

which provides that unknown parameter rgt disappears from further considerations, and 
that Ps(r) given by Eq. ( IBIOI) and Eq. ( IBllI) really satisfies the equation f l3.12p . Then, from 



Eqs. flBlOp and fIBllI) it follows the sought expression (13.141) for the member ne-s{r). Finally, 
from Eqs. (I3.1ip and (13.141) it is obtained the other condition which is necessary for the 
determination a and b, namely 

a ■ exp {-Ko-ers;i) + b ■ exp (Ko;e?^s;i) = 1 - a, (B16) 
where the free parameter a has to be determined from the condition Eq. (13.91) . 



APPENDIX C: THE SYSTEM S, 



The region: rg-e < r < oo. In the case (e) we will start from the fact that the fixed probe 
particle with the charge (— e) principally can not be represent of a free electron in the system 
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Sin from the aspect of the interaction with positive ions, since their average distribution in 
the neighborhood of such a particle would resemble to a distribution of positive ions in 
plasma in the neighborhood of a heavy negative ion. 

The exit from this situation is the treatment of accessory system Si^^ as a single com- 
ponent one with the electron gas on the corresponding positive charged background, which, 
contrary to the background described in Part 1, is not homogeneous. We keep in mind the 
background with the charge density taken in the form ZiC ■ nf\r), which in the classical 
case is able to model the average distribution of positive charge in the neighborhood of an 
electron in the system Sin- If sought distribution is found, the factor rii\r) could be treated 
as the corresponding ion density. Here we use the fact that such a distribution can be found 
since existing conditions which establish correspondence between the systems S'a^'' and Sin 
are enough for determination of all needed parameters. 

One of the mentioned conditions is that the behavior of rii\r) and nf\r) reflect the 
existence of the electron-ion correlation (discussed in Section IIII Al) which is characterized 
by the coefficient a. We will take into account that this coefficient in the case (e) can be 
treated as the probability of the following event: the decreasing of number of electrons for 
Zi in the region r^.f. < r < oo coincides with the decreasing of the number of ions for unity. 
This means that in this case we have the relation 



which connects rie (r) and nf {r) in the region r^je < r < oo by means of the known param- 
eter a. In Section nil BI the relation (ICip is taken as the starting point for the determination 
nf\r) and nt\r) outside of probe particle self-sphere in the system S^f^ . 

The region: < r < Tg-^e- In order to determine the member ni-s{r) in Eq. fl3.26p we 
will take into account the following facts: the form of nj.s(r) has to be similar to the form 
of the member ne;s(r) in the case (i); the parameters which characterize ^^.^(r) have to be 
closely connected with the parameters which characterize ne;s(r); the procedure of obtaining 
of ni-s{r) has to provide automatic applicability for all possible Ng, Zi, and Tj and self- 
consistence of final expressions. In accordance with this, the member rii-sir) will be taken 
as a superposition of two dimensionless functions of dimensionless argument (r/rs-e), given 
by relations 



a = 




(CI) 



rii-Ar) =ai- {r/r^-e) ^ exp [-(K^.ir^.e)r/r^;e] + h ■ {r/rs-e) ^ exp [(fi;^.ir^.e)r/r^;e], 



(C2) 
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^ S'i ~ 

f^s;i = f^O;e " '- = f^O;e " , (C3) 

where kq-s is given by (13.131) . and the screening constant Hg-i is chosen in such a manner 
which guarantees that the relation 

r' , r' . r" , r" , 

exp {Tfis;irs;e ) = GXp (=F/«0;e'"s;i ) (C4) 

S]i ' S]i ' s;e ' s;e 

is vahd for any (r'/rg-i) = {r" /rs-e), when < r' < Vg-i and < r" < r^^g. The coefficients 
and 6j in Eq. (1020 have to be found from the conditions (13.271) and (I3.28p . 



APPENDIX D: THE INTERPRETATION OF THE SYSTEMS CONSIDERED 

Let Sin-Mi tie the model finite spherical system with total ion and electron numbers Mj 
and Me = ZjMj, where Mj is an integer number, and with the radius and volume 
Vm, determined by relation: Mi/Vui = Ni. From here it follows that the system Sin-M, is 
neutral as a whole and Me/Vui = Ng = ZiNi. We will assume that in Sin-M^ the ion and 
electron components can be treated as gases in states of the thermodynamical equilibrium 
with temperatures Tj and Tg. In a usual way, we will treat the basic system as a 
thermodynamical limit of the systems Sin-M^ , i-e. as the result of transition: Mj oo and 
Vm, — > oo, under conditions 

Mi/VM, = Ni, Ti = const., = const. (Dl) 

In the case (i), we will associate with every system Sin- Mi an other system S'^^\j^, which 
differs from Sin-^Ui only by the change of one of the free ions for the probe particle, with the 
same charge ZiC fixed in the center of that system. 

We will take into account that the systems 5"^*]^/, are also neutral as a whole. From here 
follows the relation 
Rhii 

ZiC + Zie- j ni-^M, (r) ■ Ar^r^dr - e ■ N,Vm^ = 0, (D2) 



where ni-Mi{f) are the corresponding mean local ion densities, and 

A-jr r 

N,Vm, = Z,N, ■ —RIj^ = Z,- N,- Anr^dr. (D3) 
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From (lD2p and (lD3p . after their multiplication with (Zic) ^ the equation follows 

Km, 

[A^i-ni;M.(r)]-47rr2dr = 1. (D4) 



Based on it, we have it that 



Rm 



lim / [Ni - Ui.MXr)] ■ Anr^dr = [Ni - rii-ocir)] ■ ATrr^dr = 1, (D5) 

Mi^oo J J 


where ni.^{r) = lim rii-Miif') for any r > 0. Here, we will treat just rie-ooif^) as the mean 
local ion density in the accessory systems Sa \ i.e. as nf\r). In accordance with this, from 
fIDSP it directly follows that n[^\r) has to satisfy the condition fl2.13p for the case (i). 

Since a similar reasoning may be repeated in the case (e), we can consider the condition 
fl2.13p must be satisfied in both (i) and (e) cases. Finally, it is clear that in the case when 
one of components in described systems is changed with the corresponding background, we 
obtain the equation which corresponds to Eq. (21) from Part 1. 
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